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Progress in Three-Dimensional Coherent X-Ray Diffraction Imaging 
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The Fourier inversion of phased coherent diffraction patterns offers images without the resolution 
and depth-of-focus limitations of lens-based tomographic systems. We report on our recent experi- 
mental images inverted using recent developments in phase retrieval algorithms, and summarize ef- 
forts that led to these accomplishments. These include ab-initio reconstruction of a two-dimensional 
test pattern, infinite depth of focus image of a thick object, and its high-resolution (~ 10 nm res- 
olution) three-dimensional image. Developments on the structural imaging of low density aerogel 
samples are discussed. 
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I. INTRODUCTION 

In the last five years or so several new ideas have 
combined to provide us with a working solution to the 
phase problem for non-periodic objects. This capabil- 
ity opens exciting possibilities for using coherent x-ray 
diffraction microscopy (CXDM) for 3D imaging of few- 
micron-sized objects at resolution levels previously in- 
accessible to x-ray microscopy. Since the first proof of 
principle demonstration of CXDM^ , a number of groups 
have been working to bring these possibilities into real- 
ity. Recent estimates^ of the dose and flux requirements 
of such measurements, indicate that attractive resolution 
values (about 10 nm for life science and 2-4 nm for mate- 
rial science) should be possible with reasonable exposure 
times using modern synchrotron-radiation sources. Thus 
CXDM promises a 3D resolution limited only by radi- 
ation damage, the collection solid angle and the num- 
ber of x-rays collected. We therefore expect to have an 
advantage over lens-based tomography schemes that are 
currently limited in resolution and efficiency by the lens 
fabrication technology and, in principle, by the depth 
of focus effect. This capability provides an extremely 
valuable tool for understanding nanoscience, such as the 
study of minimum energy pathways for crack propagation 
in brittle solids, and characterizing the internal structure 
of mesoporous structures that are synthesized for a wide 
range of applications. 

In this paper we review the historical developments 
which have led to these opportunities and describe some 
of the activities of our multi-institutional collaboration, 
working at beam line 9.0.1 at the Advanced Light Source 
at the Lawrence Berkeley National Laboratory. In par- 
ticular we will describe here two experiments, which 
demonstrate spectacular 3D imaging at 10 nm resolution. 
The portion of the work devoted to life-science imaging. 



largely by the Stony Brook group, has been reported in 
these proceedings by Limai^ 



II. CONCEPTUAL HISTORY 

The observation by Sayre in 1952"' that Bragg diffrac- 
tion under-samples the diffracted intensity pattern was 
important and led to more specific proposals by the 
same author for X-ray diffractive imaging of non-periodic 
objects.— These ideas, combined with the rapid develop- 
ment of computational phase retrieval in the wider op- 
tics community especially the "support constraint"i^i2i^, 
enabled the first successful use of CXDM. An important 
review, which attempted to integrate the approaches of 
the optical and crystallographic communities, appeared 
in 1 990 »^ The connection was made between the "solvent- 
fiattening" or " density- modification" techniques of crys- 
tallography and the compact support requirements of the 
hybrid input-output (HIO) algorithm. The importance 
of fine sampling of the intensity of the measured diffrac- 
tion pattern was recognised at an early stage^^^ and has 
led to the method being referred to as " oversampling" , 
since the Shannon sampling interval is half the Bragg 
interval in each dimension. The result of Shannon sam- 
pling the intensity is that the diffracted phased amplitude 
is at least two-fold oversampled in each dimension, which 
implies that the object obtained by transformation of an 
exactly-known amplitude function will be surrounded by 
a zero-padded region of at least three times the object 
area for 2D reconstructions, or at least seven times the 
object volume for 3D reconstructions. Such zero padding 
is a necessary concomitant of the use of a support con- 
straint, although it is sufficient to oversample at less than 
the factor of two per dimensiorkii. The use of the support 
constraint as a means to phase the diffraction pattern is 



2 



now very widespread in the growing CXDM community. 
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FIG. 1: SEM images of gold ball clusters (left) and recon- 
structed soft X-ray image (right) recorded at A = 2 nm^ 



III. EXPERIMENTS WITH X-RAYS 

In spite of the promise shown in simulations, the above 
theoretical advances were not accompanied by immedi- 
ate experimental progress in the practical application of 
phase retrieval. The first successful X-ray technique was 
developed by the Stony Brook group at the XI undulator 
beam line at the National Synchrotron Light Source at 
Brookhaven. The fruit of this effort, reported by Miao, 
Charalambous, Kirz and Sayre in 1999ji was the first in- 
version of an experimental X-ray diffraction pattern to 
an image of a non-periodic object at 75 nm resolution. 
This success proved to be the beginning of a significant 
expansion in interest in CXDI in the US. 

In the last few years CXDM activities in the US has 
involved four groups which have all made contributions 
to the XRM 2005 conference: Stony Brook/Brookhaven, 
University of California at Los Angeles (UCLA) , Univer- 
sity of Illinois / Argonne, and ourselves at University of 
Arizona / Livermore Lab / Berkeley Lab. Stony Brook / 
Brookhaven have constructed a sophisticated experimen- 
tal station for tomographic imaging of life-science speci- 
mens at cryogenic tempcraturest^ This apparatus is now 
installed at ALS beam line 9.0.1— and serves all of the 
groups doing CXDM at the ALS. Robinson and cowork- 
ers at the University of Illinois have applied the princi- 
ples of CXDM to hard x-ray experiments on microcrys- 
talline particles, the density variations of which produce 
a diffraction pattern centered on each Bragg spot. The 
pattern can be reconstructed in 2D^^ or scanned in 3D by 
very slight rotation of the crystal to give the equivalent 
of a tilt series. Such data have been reconstructed tomo- 
graphically to produce a 3D image at 80 nm resolutionpi^ 
Miao (now at UCLA) and coworkers have made consid- 
erable progress in pushing the CXDI method to higher 
resolution in 2D (7 nm), higher x-ray energies and to 



FIG. 2: Scanning Electron Microscopy (SEM) of a silicon 
nitride pyramid (right) and the gold balls deposited on the 
hollow side of the membrane (left). Scalebar is 1 nm. 



a limited form of 3D>iS, They have also made the first 
application of CXDM to a biological sample»ii 

Our own efforts in this area have concentrated on re- 
trieving phase based on the diffraction data alona"^, three 
dimensional ab-initio reconstruction of a test 3D ob- 
ject made of 50 nm gold balls deposited on a pyramid 
shaped silicon nitrade membrane at 10 nm resolution^ 
two dimensional images with infinite depth of focus, and 
the imaging of materials sciences samples such as Ta205 
aerogel foams. 



IV. EXPERIMENTAL RESULTS 

Our experiments reported here were carried out us- 
ing the abovementioned Stony Brook University diffrac- 
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FIG. 3: (a) Reconstruction from a single view diffraction pat- 
tern at an object orientation ip = 24°. Scalebar is 1 ^m. (b) 
Infinite depth of focus projection images, for the object orien- 
tation (fi — 0° . (c) 2D section of the three dimensional diffrac- 
tion pattern: as the object is rotated, the recorded Ewald 
spheres intercept the plane leaving some gaps of missing data. 



3 




tion apparatus, as well as an earlier apparatusaS. In the 
Stony Brook apparatus we carried out experiments using 
750 eV (1.65 nm wavelength) X-rays that were selected 
from the undulator spectrum by a zone-plate monochro- 
mator with a spectral resolution of A/AA=:1000. The 
4-//m-diameter monochromator exit pinhole also selects 
a transversely spatially coherent patch of the beam. The 
sample was located 20 mm from this pinhole. A direct- 
detection bare CCD detector, with 20 /ini pixel spac- 
ing, 1340x1300 pixels, was located 142 mm behind the 
sample. At these CCD and wavelength settings we have 
an object sampling interval in x and y of Ax = 9.8 nm 
(in the small-angle approximation) and a field width of 
w — NAx= 11.7 /im. A beam-stop blocks the direct 
undiffracted beam from impinging on the CCD. More 
details are given by Beetz et alJ^^. Diffraction patterns 
were collected with the sample oriented at rotation an- 
gles (f> of -60° to -1-60°, at 1° intervals. Total exposure 
time was about 3 hours per sample. 

Test samples were made by placing a droplet of so- 
lution containing unconjugated colloidal gold balls on a 
silicon nitride membrane (thickness 100 nm) and allowing 
it to dry. A two dimensional object was imaged without 
prior knowledge about its shape by periodically updating 
the support region based on the current object estimate 
(Fig. 

A three-dimensional test sample was produced by 
placing a droplet of colloidal gold solution on a three- 
dimensional silicon nitride pyramid-shaped membrane^ 
This drop quickly evaporated and left the gold balls in 
a characteristic pattern where the gold tended to fill in 
the edges of the pyramid. An SEM image of the object is 
shown in Fig. [21 The pyramid base width is 2.5 /xm and 
the height (base to apex) is 1.8 //m. An earlier, larger, 
silicon nitride pyramid object is shown on the right side 
of Fig. El 

Two dimensional projection images may be recovered 
from the diffraction intensities without having to first 
undergo a full 3D reconstruction, and we found this is 
a useful step to quickly examine our 3D datasets. The 
diffraction intensities from a single sample orientation are 
recorded on the Ewald sphere and will have the same 
depth of focus as a microscope with NA equivalent to the 
solid angle intercepted by the CCD. For our experimental 
parameters, giving NA = 0.084, we have a depth of focus 
of 120 nm, which is considerably smaller than the 1.8 
/im thickness of the pyramid object. A reconstructed 
image, from a single- view diffraction pattern is shown in 
Fig. jSl^a). Artifacts due to defocus are clearly seen in 
the image. By the choice of the parabolic term of the 
retreived phases^^ the plane of focus can be numerically 
scanned throughout the depth of the object. 

By sectioning the three dimensional Fourier space with 
a sphere of larger radius than the Ewald surface, we can 
increase the depth of focus. Infinite depth-of- focus two- 
dimensional projection images were obtained (Fig ^p) 
from a plane in reciprocal space perpendicular to the pro- 
jection direction (Fig. |31;)i2i. 




FIG. 4: Three dimensional diflraction pattern (left) (with 
a quadrant removed for visualization) and reconstructed 3D 
images— (right) showing the isosurface as well as the projec- 
tion images of the sample. 




FIG. 5: (a) SEM of an areogel sample with reference points 
used to help the phase retrieval, and (b) two dimensional 
projection of the reconstructed 3D image at 14 nm resolution. 

A full 3D image was obtained^^ by performing phase 
retrieval^i^i on the entire 3D diffraction dataset. The 
reconstructed volume image reveals the structure of the 
object in all three dimensions and can be visualized in 
many ways including isosurface renderings, projections 
through the data (Fig. or slices (tomographs) of the 
data. 

We have applied 3D diffraction imaging to determin- 
ing the 3D structure of low density aerogel foam samples. 
These metal oxide foams are low density (100 mg/cc) and 
have an internal skeleton structure composed of Ta2 05. 
Our full 3D reconstructions were obtained with no a pri- 
ori information about the sample, including no measure- 
ment of the missing low spatial frequency data in the 
beamstop region. The reconstructed image, shown in 
Fig. 13 reveals not only the particle shape, but also in- 
ternal foam structure such as the strut geometry, which 
can be used to calculate the foam mechanical properties. 
Further details will be given in another paper. 



V. PHASE RETRIEVAL AND IMAGE 
ANALYSIS 

We overcame two key computational challenges in im- 
plementing high-resolution 3D phase retrievaljSS specif- 
ically performing the numerous Ik'^ FFTs required for 
phase retrieval in a reasonable time and managing the 
memory requirements of large 3D data sets. Memory 
and calculation requirements are significant and suggest 
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a cluster-based solution. The FFTs (distJFt) have been 
optimised for the G5 vector processor architecture by 
the Apple Advanced Computation Group and uses stan- 
dard MPI interfaces to perform distributed giga-elcmcnt 
or larger FFTs. Reconstruction code is written in C, is 
fully parallelised, and uses distributed memory and MPI 
interfaces to share the workload across all CPUs in the 
system. This includes application of real and Fourier 
space constraints and dynamic support refinement using 
the Shrinkwap algorithm. 

Using 16 G5 Macintosh computers with Infiniband in- 
terconnects we achieve an FFT speed of 7.9 sec/FFT on 
a 1024^ voxel data set, giving a full reconstruction time 
of 14 hours (based on 2000 iterations, 2 FFTs per itera- 
tion plus other floating point operations needed for the 
reconstruction). Timings for a 512'^ data cube are 850 
rosec/FFT, enabling us to perform a full reconstruction 
in 1.5 hrs. 



agreement of the Fourier amplitudes of the reconstructed 
image with the data, back up this assertioniSii^. 

VI. HOLOGRAPHIC-ENHANCED PHASE 
RETRIEVAL 

It was noted that the autocorrelation functions in 
some of this work also included faithful, although low- 
resolution, holographic images of some of the clusters, 
due to the occurrence of a single isolated ball near the 
object. In analogy with the "heavy atom" method of 
crystallography, by placing a reference point object near 
the sample we can obtain a one-step, estimate of the sup- 
port function. Although the holographic image is noisier 
than the recovered image (Fig. 0, it provides a useful 
starting point to the algorithm. 
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FIG. 6: Line-outs of the real part of the reconstructed com- 
plex amplitude 3D image, through the isolated single ball at 
the pyramid apex. Dashed lines show a simulated 3D coherent 
image with a cube OTP of 10 nm resolution and with a 60° 
missing sectoi^. The lineout through three lines (d) demon- 
strates that the resolution is sufficient to clearly distinguish 
between different 50 nm gold spheres. 




FIG. 7: Fourier transform of the diffraction pattern of the 
object illuminated with a large beam: the central part of the 
picture shows the autocorrelation of the pyramid, but some 
reference points produce off-centered holograms. 

Inspired by this holographic method to help the phase 
retrieval step (see also^^, we developed a methodology 
to deposit controlled reference points near the object by 
metallorganic deposition using a focused ion beam. Our 
initial trials arc illustrated in Fig. |S1 
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Although we cannot exactly quantify the resolution of 
the image, which would require knowing the object's 3D 
structure, our analysis shows we can consistently retrieve 
phases out to the maximum spatial frequency recordediS 
(further improvements in.'^^ A line-out through the re- 
constructed image can easily resolve 50 nm spheres that 
are touching each other (see Fig. EJ. From such image 
line-outs, and comparisons of reconstructed X-ray im- 
ages with the SEM image, we have confidence that our 
achieved image resolution is close to 10 nm. Further anal- 
ysis of the consistency of the retrieved phases, and the 




coccolith 



FIG. 8: SEM image of a coccolith shell deposited on a silicon 
nitride membrane (left). A reference point deposited with a 
focused ion beam produces a holographic image (right). 
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VII. CONCLUSIONS 

We have demonstrated ab-initio two dimensional im- 
ages with infinite depth of focus, and three dimensional 
images of test objects at a resolution of 10 nm. Pre- 
liminary images of areogel foams were presented. These 
images of complicated and unknown objects, along with 
the rigorous analysis of known test objects, show the ro- 
bustness of our ab initio phase retrieval technique. In the 
case of the aerogel particle, the reconstruction was per- 
formed "blind" without the operator (A. Barty) of the 
reconstruction software aware of the SEM image, or the 
size or shape of the object. 

While the recent experimental progress to date has 
been rapid and extremely encouraging, we are looking 
forward to further improvements in the technique, includ- 
ing faster acquisition times (with an improved beamline) 
that will allow us to achieve even higher image resolu- 
tion. Given the scaling of required dose to the inverse 
fourth power of resolution^, and estimates of coherent 
flux improvements achievable with an optimized beam- 
line and undulator source, we estimate that we should be 
able to achieve resolutions of 2-4 nm on material science 
samples. Possible^ applications for the technique include 
characterizing the pore structure of vesicular basalt, the 
formation of voids in metals, and many other investiga- 



tions of the nanoworld. The techniques that we have 
developed will also be applied to exciting new prospects 
for imaging of large macromolecules and assemblies at 
near atomic-resolution imaging, which will be achieved 
using X-ray free-electron lasers^i and aligned molecule 
diffractionSi. 
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